Astrophysical constraints on the confining models : the Field Correlator Method 
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We explore the relevance of confinement in quark matter models for the possible quark core of 
neutron stars. For the quark phase, we adopt the equation of state (EoS) derived with the Field 
Correlator Method, extended to the zero temperature limit. For the hadronic phase, we use the 
microscopic Brueckner-Hartree- Fock many-body theory. We find that the currently adopted value 
of the gluon condensate G2 — 0.006 — 0.007 GeV , which gives a critical temperature Tc — 170 MeV, 
produces maximum masses which are only marginally consistent with the observational limit, while 
larger masses are possible if the gluon condensate is increased. 
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nuclear matter one is the main point to understand if 
a core of pure quark matter can exist in NS. This pos- 
sibility has been addressed and extensively discussed in 
the literature [a, [a, 0]- I^ the NJL model, where the 
phase transition corresponds to chiral symmetry restora- 
tion, the quark onset at the center of the NS, as the 
mass increases, marks an instability of the star, i.e. the 
NS collapses to a black hole at the transition point since 
the quark EoS is unable to sustain the increasing central 
pressure due to gravity. Indeed, at the maximum mass 
the mass-radius relation is characterized by a cusp [g]. 
On the contrary, for the quark EoS based on the MIT 
bag model, it is possible to find a range of the various 
parameters which corresponds to a stable NS. It must 
be noted that stability is also found in other approaches 
that explicitly take into account the dynamics of confine- 
ment, such as the dielectric model 9] or a modification 
of the NJL with an ad hoc confining potential [i3|- A 
modified NJL with the explicit inclusion of colour su- 
perconductivity and isoscalar vector meson coupling |ll| 
produces stable NS as well. This shows how the pres- 
ence of quark matter in the interior of NS depends on 
the adopted quark matter model. 

The intriguing relation between the stability of NS 
with a quark matter core and confinement has already 
been addressed in [12| , and in this paper we elaborate fur- 
ther on this idea by resorting to the EoS of the quark mat- 
ter at finite temperature and density, obtained in the non 
perturbative framework of the Field Correlator Method 
(for a review see [4|), which gives a natural explanation 
and treatment of the dynamics of confinement in terms 
of Color Electric (CE) and Color Magnetic (CM) field 
correlators. In this way the FCM will be tested by com- 
paring the results for the neutron star masses with the 
existing phenomenology, which turns out to be a strong 
constraint on the parameters used in the model. 

It will be shown that this approach, unlike the non- 
confining NJL model, admits stable NS with gravita- 
tional masses slightly larger than 1.44 M©, and this ap- 
plication of the FCM to the study of NS, which has not 
been considered before, provides definite numerical indi- 
cations on some relevant physical quantities, as the gluon 



I. INTRODUCTION 

QCD at finite temperature and density is the essen- 
tial theoretical tool to describe various interesting phe- 
nomenological sectors from relativistic heavy ion colli- 
sions to the inner structure of neutron stars. In the large 
temperature and small density region both experiments 
[i| and lattice simulations [2| clearly indicate that the 
theory is in a non-perturbative regime at least up to 
temperatures T ~ 3rc(0) (Tc{0) is the deconfinement 
temperature at zero quark chemical potential iiq = 0). 
However, in the opposite region of the phase diagram, 
i.e. at small T and large fig, where strong coupHng ef- 
fects are expected as well, no QCD lattice simulations 
are available yet. 

Due to the lack of lattice data, analytic approaches 
based on more elementary models, such as the Nambu- 
Jona-Lasinio (NJL) model Q, that mimic some non- 
perturbative features of QCD, are mostly used in the 
large density region, typical of neutron star interior. Un- 
fortunately the NJL model cannot be used in the other 
limit of zero chemical potential and high temperature be- 
cause of the lack of the gluon degrees of freedom. This is a 
general feature of many models, which cannot make pre- 
dictions for both limits, i.e. high temperature and zero 
chemical potential or high chemical potential and low 
temperature. This is clearly a serious drawback, since 
the models cannot be fully tested. One of the few excep- 
tions is the Field Correlator Method (FCM) Q, which in 
principle is able to cover the full temperature-chemical 
potential plane. Furthermore the method contains ab 
initio the property of confinement, which is expected to 
play a role, at variance with other models like, e.g., the 
NJL model. 

The study of the properties of neutron stars (NS) con- 
cerns the large density (and low temperature) region of 
the phase diagram and, in particular, it requires the QCD 
non-pertubative Equation of State at small T and large 
fiq. The comparison of the quark matter EoS with the 



condensate, to be compared to the ones extracted from 
the determination of the critical temperature of the dc- 
confinement phase transition. This shows the relevance 
of the comparison of the model predictions in the high 
chemical potential region with the astrophysical phe- 
nomenology, which is one of the main purposes of the 
present paper. 

In the next section the FCM at finite temperature and 
density is briefly recalled, while Sec. Ill contains some 
details of the EoS for the hadronic phase. Our analysis 
of stability of the NS is presented in Sec. IV and finally 
Sec.V is devoted to the conclusions. 



II. QUARK MATTER: EOS IN THE FIELD 
CORRELATOR METHOD 



A systematic method to treat non perturbative effects 
in QCD is by gauge invariant field correlators [j|. The 
approach based on the FCM provides a natural treatment 
of the dynamics of confinement (and of the dcconfinement 
transition) in terms of the CE {D^ and Df) and CM 
{D^ and D^) Gaussian (i.e. quadratic in the tensor 
F""^) correlators. D'-^ and Df are related to the simplest 
non trivial 2-point correlators for the CE and CM fields 
by: 

g^ < Trf[a{xMx,y)Ck{yMy,x)] >= 

5,k[D^ {z) + D\ {z) + zl ^; ' ]±z,zk ^; (1) 



dz 



dz^ 



where z — x — y, C indicates the CE (E) field or CM 
(H) field (the minus sign in the previous expression cor- 
responds to the magnetic case) and finally 



$ = P exp 



W 



A^'dz,, 



(2) 



is the parallel transporter. 

D'~^ and I?f have a perturbative contribution which is 
responsible of their singular behavior at z ~ {D ~ z^** 
for z — > 0), but also a non-perturbative part which is 
normalized to the gluon condensate [1|. D^ contributes 
to the standard string tension and is directly related to 
confinement, so that its vanishing above the critical tem- 
perature implies deconfinement. 

The FCM has been extended to finite temperature T 
and chemical potential fiq in order to describe the de- 
confinement phase transition IJ, Hi Ej Ea Ei Il3l ■ ^^ 
particular, at fig = the analytical results in the gaus- 
sian approximation, valid for small vacuum correlation 
leng ths, are in reasonable agreement with lattice data 
[11 [H, [3. The extension in ref. [H of the FCM to 
finite values of the chemical potential, allows to obtain a 
simple expression of the Equation of State of the quark- 
gluon matter in the relevant range of baryon density. The 
comparison of this EOS with a realistic baryonic EOS will 
be the crucial point of our investigation. 



It must be noticed that the generalization of the FCM 
at finite T and /ig provides an expression of the pres- 
sure of quarks and gluons where the leading contribu- 
tion is given by the interaction of single quark and gluon 
line with the vacuum, called Single Line Approxima- 
tion (SLA), while the pair and triple correlations yield 
higher order corrections. In the SLA, within few per- 
cent, the auark pressure, for a single flavour, is given by 

MMMM 



Pg/T =-^[M 7^ 

where v — niq/T ^ and 
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and Vi is the large distance static QQ potential: 

r-l/T 



Fi=/ dril-TT) dxX^f(VF+^) (5) 
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The gluon contribution to the pressure is 

1 



''«/^' = s? 



c^XX' 



exp(x + |^)-l 



(6) 



Note that the potential Vi in Eq.(l5]) does not depend 
on the chemical potential and this is partially su ppo rted 
by lattice simulations at small chemical potential [la.[l9|. 
In our opinion, although we are considering the range 
T ^ (in the following calculations we fix the value 
T = 1 MeV) and large Hq, relevant for the NS, this ap- 
proximation is still valid. Indeed the non perturbative 
contribution to Df{x) is parametrized as 



i?f(x)=i?f(0)exp(-N/A) 



(7) 



where A is the correlation length (0.34 fm for full QCD) 
and the normalization is fixed by the condition at T = 



D^iO) + Df{0) = -G, 



(8) 



G2 is the gluon condensate whose numerical value, de- 
termined by the QCD sum rules, is known with large 
uncertainty [2(j | 



G2 = 0.012 ± 0.006 GeV 



(9) 



According to [16| , the critical temperature at /i = in 
the FCM turns out to be T ~ 170 MeV for G2 -- 0.006 
GeV^. If confinement is dominated by non-perturbative 
contributions, the normalization Df (0) in Eq. ([7]) can 
be indeed identified with the term appearing in Eq. ([5]), 
which has been denoted by the same symbol. Then from 
Eqs. dni, 0, dH), in the limit T ^ 0, we get 



Vi{T^Q) < ^GzA^ 



TT 
"9 



(10) 



However other choices of Vi are possible, and these will 
be considered in the Discussion section. 

Since, on general ground, we expect that the value of 
the gluon condensate decreases at large densities 21], the 
assumption that Vi is /i- independent, should not quali- 
tatively modify our analysis. 



III. HADRONIC PHASE: EOS IN THE 
BRUECKNER-BETHE-GOLDSTONE THEORY 

The EOS constructed for the hadronic phase at 
T = is based on the non-relativistic Brueckner-Bethe- 
Goldstone (BBG) many-body theory [22, which is a 
linked cluster expansion of the energy per nucleon of nu- 
clear matter, well convergent and accurate enough in the 
density range relevant for neutron stars. In this approach 
the essential ingredient is the two-body scattering matrix 
G, which, along with the single-particle potential f/, sat- 
isfies the self-consistent equations 

U{k-p) = J2 {kk'\Gip;e{k) + e{k'))\kk')a (12) 

where v is the bare nucleon-nuclcon (NN) interaction, p 
is the nucleon number density, to is the starting energy, 
and \kakb)Q{kakt\ is the Pauli operator. e(fc) = e(fc; p) — 
■^k^ + U{k; p) is the single particle energy, and the sub- 
script "a" indicates antisymmetrization of the matrix ele- 
ment. In the BHF approximation the energy per nucleon 
is 



E , , ih^ kl 
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(13) 
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For the two-body interaction v, we choose the Argonne 
1)18 nucleon- nucleon potential [23| . We have also intro- 
duced three-body forces (TBF) among nucleons, adopt- 
ing the phenomenological Urbana model [2j|. This al- 
lows to reproduce correctly the nuclear matter saturation 
point po ~ 0.17 fm^^, E/A « —16 MeV, and gives values 
of incompressibility and symmetry energy at saturation 
compatible with those extracted from phenomenology 
[25!|. Moreover, the BBG approach has been extended to 
the hyperonic sector in a fully self-consistent way [26|, [231 , 
by including the E^ and A hyperons. 

In this paper, we adopt a conventional description of 
stellar matter, as composed by neutrons, protons, and 
leptons in beta equilibrium [2g|. The EoS for the beta 
equilibrated matter can be obtained for a given composi- 
tion, together with the chemical potentials of all species 
as a function of the total baryon density. The chemical 
potentials are the fundamental input for the equations of 



chemical equilibrium, charge neutrality conditions, and 
baryon number conservation, i.e.. 



Mn = 



Me 



Pp + Pe- 

Mm- 



Xr, 



x^ 



p 



(15) 
(16) 

(17) 
(18) 



where Xi — pi/ p is the nucleonic fraction of the species 
i. The above conditions allow the unique solution of a 
closed system of equations, yielding the equilibrium frac- 
tions of the nucleonic and leptonic species for each fixed 
nucleon density. Once the composition of the /3-stable, 
charge neutral stellar matter is known, one can calculate 
the equation of state, i.e., the relation between pressure 
P and energy density e as a function of the baryon density 
p. It can be easily obtained from the thermodynamical 
relation 



P 
Pb 



dE 



= P 



dV 

2d{eB/p) 



= Pb + Pi 



dp 



Pi=p- 



djei/p) 
dp 



(19) 
(20) 



with E the total energy and V the total volume. The 
total nucleonic energy density is obtained by adding the 
energy densities of each species e^. As far as leptons are 
concerned, at those high densities electrons are a free ul- 
trarelativistic gas, whereas muons are relativistic. Hence 
their energy densities e/ are well known from textbooks 

M. 



IV. PHASE TRANSITION IN BETA-STABLE 
MATTER AND NEUTRON STAR STRUCTURE 

We are now able to compare the pressure of the two 
phases, namely the pressure in the hadronic phase given 
in Eqs. p9)) . ((20)) with the one in the quark-gluon phase 
which, according to [la, 13' '^^^ t)e written as 



Pgg ^ Pg 



E 

j—u,d,s 



pj 

q 



Ae 



(21) 



where Pg and P-j are respectively given in Eq. ^ and 
©, and 



(11-|JV/) G2 
32 2 



(22) 



corresponds to the difference of the vacuum energy den- 
sity in the two phases, being Nf the fiavour number. 

By assuming a first order hadron-quark phase tran- 
sition [30| in beta-stable matter, we adopt the simple 
Maxwell construction. The more general Gibbs construc- 
tion [5|] is still affected by many theoretical uncertainties 
[3lj, and in any case the final mass-radius relation of 
massive neutron stars [32] is slightly affected. 
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FIG. 1: Pressure as a function of the baryon chemical poten- 
tial. The full line represents the BHF calculations, and the 
dashed ones the model discussed in this paper with two dif- 
ferent choices of the parameter Vi, and several values of the 
gluon condensate G2. See text for details. 



We impose thermal, chemical, and mechanical equi- 
librium between the two phases. This implies that the 
phase coexistence is determined by a crossing point in 
the pressure vs. chemical potential plot, as shown in 
Fig. [H There we display the pressure P as function of 
the baryon chemical potential fiB for baryonic and quark 
matter phases. In the upper panel we show the results 
obtained using Vi — 0, whereas in the lower panel cal- 
culations with Vi = 0.01 GeV (according to the indica- 
tion of the constraint in Eq. (fTO|) ) are displayed. The 
solid line represents the calculations performed with the 
BBG method with nucleons, and the other lines repre- 
sent results obtained with different choices of the gluon 
condensate G2. We recall that the chosen values of G2 
give values of the critical temperature in a range between 
160 and 190 MeV 0. 

We observe that the crossing point is significantly af- 
fected by the value of the gluon condensate, and only 
slightly by the chosen value of the potential Vi . More- 
over, with increasing G2, the onset of the phase transition 
is shifted to larger chemical potentials. Hence, we expect 
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FIG. 2; Pressure as a function of the baryon density, nor- 
malized with respect to the nuclear matter saturation density 
Po. 



that the neutron star will possess a thicker hadronic man- 
tle with increasing G2. 

In Fig. [2] we display the total EoS, i.e. the pressure 
as a function of the baryon density for the several cases 
discussed above. The plateaus are consequence of the 
Maxwell construction. Below the plateau, /3-stable and 
charge neutral stellar matter is in the purely hadronic 
phase, whereas for density above the ones characterizing 
the plateau, the system is in the pure quark phase. 

The EoS is the fundamental input for solving the well- 
known hydrostatic equilibrium equations of Tolman, Op- 
penheimer, and Volkov [23| for the pressure P and the 
enclosed mass m 
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being e the total energy density (G is the gravitational 
constant). For a chosen central value of the energy den- 
sity, the numerical integration of Eqs. (P5)) and (|M|) pro- 
vides the mass-radius relation. For the description of 
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FIG. 3: The gravitational mass (in units of the solar mass) is displayed as a function of the central baryon density, normalized 
with respect to the nuclear matter saturation density po (left panel), and the corresponding radius (right panel). 
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FIG. 4: Same as Fig.3 for Vi = 0.01 GeV. 



the neutron star crust, we have joined the equations of 
state above described with the ones by Negele & Vau- 
therin ?)^ in the medium-density regime (0.001 fm~ < 
p < 0.08 fm~'^), and the ones by Feynman, MetropoUs, 
& Teller ^ and Baym, Pethick, & Sutherland 35] for 
the outer crust (p < 0.001 fm~ ). 

In FigO we display in the left panel the gravitational 
mass (in units of solar mass Mq = 2 x 10'^^ g) as a function 
of the central baryon density (normalized with respect to 
the saturation value), and the corresponding radius in the 
right panel. We observe that the value of the maximum 
mass spans over a range between 1.4 and 1.8 solar masses, 
depending on the value of the gluon condensate G2, as 
shown in Table 1. The stability of the pure quark phase 
appears only for small values of G2, which are hardly 
in agreement with observational data. In fact, we recall 
that any "good" equation of state must give for the max- 
imum mass at least 1.44 solar mass, the best measured 
value so far [3^|. By increasing the value of G2, the max- 
imum mass increases as well, up to about 1.8 solar mass, 
but the stability of the pure quark phase is lost, and the 



maximum mass contains in its interior at most a mixed 
quark- hadron phase. By switching on the potential Vi, 
as displayed in FigHl we observe a trend similar to the 
case Vi — 0. Therefore, generally speaking we can con- 
clude that this model gives values of the maximum mass 
in any case below two solar mass, in agreement with the 
current observational data. However, the observational 
data indicate that NS with a mass of at least 1.6 solar 
masses do exist (37[ , and this puts a serious constraint on 
the value of the gluon condensate, which is not easy to 
reconcile with the value 0.006 GeV'', extracted from the 
comparison with the lattice data on the critical tempera- 
ture. This result emphasize the relevance of astrophysical 
data in testing different quark matter models. 



V. DISCUSSION AND CONCLUSIONS 

The problem of the appearance of quark matter in the 
NS core has been discussed by considering the micro- 
scopic EoS in the FCM where the dynamics of confine- 
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TABLE I: Properties of the maximum mass configuration for 
different values of the model parameters. 



ment is assumed to be a long range phenomenon. The 
results confirm the idea that confinement plays an impor- 
tant role to obtain a stable system under the gravitational 
pressure. However, in this case, pure quark matter can 
appear only for a certain range of the gluon condensate, 
which is mainly a parameter of the model. The com- 
parison with phcnomcnological data on NS masses gives 
strong constraints on the values of this parameter, which 
unfortunately are only marginally compatible with the 
range extracted by comparing the model with lattice data 
at zero chemical potential. However in this case the value 
of the large distance static QQ potential Vi turns out to 
be very small. Other choices are possible if Eq. ([7]) is as- 
sumed to be valid only at long range, while Eq. ([5]) is a 
true short range relationship. In this case the parameters 
Df(fS) in the two equations cannot be identified and may 
correspond to two different numerical values, and there- 
fore the value of Vi must be considered an independent 
parameter. In the comparison with lattice calculations 
[15J one finds a value Vi ~ 0.5 GeV at the critical tem- 
perature and for /i = 0. Besides that, the assumption of 
the independence of Vi on jjl can be questionable, it ap- 
pears in any case that the value of this parameter at high 
density and low temperature is quite uncertain. We have 
therefore varied the strength of Vi from the small values 
previously considered up to 0.5 GeV. The results for the 
EOS is reported in Fig. [5] for different values of Vi. One 
can see that the hadron-quark phase transition is shifted 
to higher values of the chemical potentials and therefore 
of the density. This can be expected just by inspection of 
the formula for the pressure, which is clearly a decreas- 
ing function of Vi. Actually already for Vi = 100 MeV 
the phase transition cannot occur in NS, which is then 
composed of baryon matter only, with a maximum mass 
around 2 solar masses. For higher values of Vi the tran- 
sition can possibly occur only at exceedingly high values 
of the density, and therefore the quark phase is irrelevant 
for NS physics. 

These results indicate once more a direct link between 
the NS quark content and the properties of deconfine- 
ment in the hadron-quark phase transition. More quan- 
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FIG. 5: Pressure as a function of the baryon chemical poten- 
tial, for different values of Vi and 6*2. The uppermost curve 
represents the BHF EoS with hyperons. 



titatively, if one considers that the well established values 
of NS masses never exceed « 1.6 solar masses, then these 
observational data constrain Vi to small values and in 
a narrow range, well below 100 MeV, in sharp contrast 
with values around 0.5 GeV extracted from lattice cal- 
culations. Despite the FCM is in good agreement with 
full QCD lattice data and is a well defined theoretical 
approach where confinement is, ab initio, the crucial dy- 
namical aspect, some refinements seem to be needed once 
the astrophysical data are considered. 

A relevant point to be clarified is the possible presence 
of hyperons, whose onset is expected to be around 2-3 
times the saturation density. In Fig 15] we plot two curves 
corresponding to the BHF EoS with and without the in- 
clusion of hyperons. As displayed, both curves coincide 
at small values of the baryon chemical potential and, af- 
ter the onset of the hyperons, the former curve grows 
faster becoming the uppermost one in the figure. 

Therefore, only at small values of Vi, of the order of 
0.01 GeV or below, the transition to quark matter oc- 
curs at about the same density as the hyperon onset, as 
displayed in Fig. [2] and [5] On the other hand, Fig[5] 
clearly shows that, when increasing Vi up to Vi ~ 0.5 
GeV, no crossing with the quark matter EOS is possible. 
We remark that the baryonic EoS with hyperons in the 
BHF framework produces a maximum mass close to 1.3 
Mq, below the observational limit, and therefore it is not 
acceptable j38| . 

It has to be pointed out that in all cases where no phase 
transition to quark matter is possible, with or without hy- 
perons, nuclear matter can reach densities where baryons 
are so closely packed that keeping their identity is highly 
questionable. This is the main physical qualitative argu- 
ment that suggests as likely a transition to quark matter. 

Another approximation used in the FCM is the so 
called Single Line Approximation where the relevant dy- 
namics is related to the interaction of a single quark or 
gluon with the vacuum. At large density this could be 



no longer true, but in the FCM the most important non 
perturbative effects are included in the field correlators 
and in particular in the gluon condensate which drives 
the transition. Therefore at large density the main ef- 
fect should be related to the fi dependence of G2 in Eq. 
([HI). Lattice data at large temperature and small den- 
sity show that the Color Electric condensate goes to zero 
at the transition point and the Color Magnetic conden- 
sate survives at large temperature. In our analysis the 
same behavior has been assumed at small temperature 
and large density (see Eqs. (pij) and (f^ ). Of course our 
results depend on this assumption and wc checked if a 



different qualitative conclusion is reached with a density 
dependent gluon condensate. Following the suggestion 
in [2l| for the density dependence of G2 we obtain the 
same qualitative results with the possibility of a larger 
NS mass ~ 2 solar masses. 
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